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A trivariate extreme value distribution 
has been derived from ihe logistic 
model for the multivariate extreme 
value distributioii. The construction of 
Its corresponding probability distribu- 
tion and density fiinctlon is described. 
In order to obtain the parameters of 
such a trivariate distribution, a general- 
ized maximum likelihood estimation 
procedure is described to allow for ihc 
cases of samples with different record 
lengths. Furthermore the reliability of 
[he estimated parameters of the irivari- 
ate extreme value distrlbuiion \s mea- 
sured ihrough ihe use of relative 



information ratios, A region in North- 
ern Mexico with six gauging stations 
has been selected to apply the trivari- 
ate model. Results produced by the 
proposed model have been compared 
with those obtained by general extreme 
value (GEV) distribution functions. 
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1. Introduction 

Flood frequency analysis has been carried out by 
using univariate distribution functions, the extreme 
value distributions being an important set of distri- 
butions used in this field of study. Generally, 
parameters of such distributions are estimated 
from a short record of flows. The variability of 
these estimates has prompted exploration of joint 
estimation models which use information from 
streamflow records of neighboring gauging stations. 

In pioneering papers Finkelstein []], Tiago de 
Oliveira [2], and Gumbel [3] gave the foundations 
for the multivariate approach to extreme value dis- 
tributions. Following this work^ several bivariate 
extreme value models began to appear in the liter- 
ature* Rueda [4] explored the logistic and mixed 
models for bivariate extreme value distributions 
when both marginals are extreme value type I 
(EVI) distributions. He reported improvements in 



the estimation of parameters when the bivariate 
approach is used. Raynal [5] developed and ap- 
plied three bivariate options from the logistic 
model of bivariate extreme value distribution for 
flood frequency analysis. He found that there exists 
an improvement in the parameter estimation 
phase^ even in the case when both samples have 
the same record lengths. 

Herein* the trivariate approach of multivariate 
extreme value distribution is presented with a view 
to its application to flood frequency analysis. 

General characteristics, the procedure for 
estimation, and reliability of parameters of the 
trivariate extreme value distributions will be de- 
scribed in the following sections. An actual applica- 
tion of the proposed model to six gauging stations 
in Northern Mexico is presented in the paper. 
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2. Characteristics of the Trivariate Logis* 
ttc Model 

From the multivariate extension of the logistic 
model for bivariate extreme value distribution [3], 
the trivariate approach is: 

F(;t;y,r,e)-exp{-[(lnF(x)r 

+ (-]n F(y)r-f (-In FCz))"}*^" , (1) 

where m is the association parameter (m^l) and 
F{s)=F(s, 0) is the marginal distribution function 
of s. Equation (1) must satisfy the following in- 
equalities (Tiago de Ohveira [6, 7]); 

F(x) F(y)F(z) < F{x, y, z) ^ min[F(x), F(y), F{z)\ 

(2) 
{F(x,y,)F{x.z)F(y,z)\^ 

-riA;y, z>_ [f^^)p^)f^^)-Y t^J 

Marginals in Eq. (1) can be either EVI distribu- 
tions: 

F(0 = exp(-exp-(^)) (4) 

or GEV distributions: 

7^(5) = exp-(l-(^))3)"^. (5) 



The combinations have been named (Escalanie 
[8]): 

a) Trivariate extreme value distribution type 111 
(TEVlll) or TriGumbel distribution. All mar- 
ginals are EVI distributions. 

b) Trivariate extreme value distribution type 112 
(TEVn2) or BiGumbel-GEV distribution, 

c) Trivariate extreme value distribution type 122 
(TEV122) or BiGEV-EVI distribution. 

d) Trivariate extreme value distribution type 222 
(TEV222) or TriGumbel distribution. All mar- 
ginals are GEV distributions. 

The particular form of Eq. (1), when the marginals 
are GEV distributions for the maxima, is (Escalante 

F{x, y, z, U], a\, j9i, u?, oi, ft, uj, tn, ft, m,) = 



b-m^) 



tfilfpt 



i^-C^hHl- (« 



where w,, a, and /S , i = 1, 2, 3, are the location, scale 
and shape parameters of the marginal GEV distri- 
butions for the maxima. The corresponding proba- 
bility density function is (Escalante [8]): 



fix, y. z, uu ctu Hu «.. en, ft. u^, m, ft. m.) =^jy^ (l" (^) ^' )""""" ' (^ " (^) ^'Y'"'" 



((l_(lZi.),.)...,(,_(Z^),,)...(,_(l^),,)..)-]. 



(7) 
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3- Estimation of Parameters 

The method of maximum likelihood for estimat- 
ing the parameters of trivariate extreme value dis- 
tributions has been chosen due to its characteristics 
for consistency in large sample estimation and ap- 
plicability in estimating the parameters of cumber- 
some density functions. 

For the case of trivariate distribution functions, 
the sample arrangements could allow having either 
an equal or different record length in any of the 
samples to be analysed. 

In order to consider all possible combinations of 
data, it is required to have a sufficiently flexible for- 
mulation, therefore the following general form of 
the likelihood function will be used based on the 
generalization obtained by Anderson [9]: 

r "3 V3 r n, -1/4 r«, I/S 



where: 



(8) 



nx.tij =are respectively the univariate and 
bivariate record lengths before the 
common period n^, 

ftj, rt5 =are respectively the bivariate and uni- 
variate record lengths after the com- 
mon period /i3. 



p =is the variable with univariate record 

before the common period, 
(p, q) -are the variables with bivariate record 

before the common period, 
(ji; y, z) = are the variables with trivariate record 

during the common period, 
(r, s) =are the variables with bivariate record 

after the common period^ 
r =is the variable with univariate record 

after the common period, 
ii =are indicator numbers such that: 

/; = 1 if /I, >0 and A=0 if n,=0. 

The logarithmic function will be used instead of 
the likelihood function. So, Eq. (8) is transformed 
into: 

LL{x,y,z, 0=/.[2;in/(pM^i)] 

+ /4[2 ln/(r., J,, eo] + /5[i ln/(r;, &)]. 

(9) 

The maximum likelihood estimators of parame- 
ters for the trivariate extreme value distributions 
are those values for which Eq. (9) is maximized. 

The corresponding logarithmic likelihood func- 
tion for the trigeneral extreme value (TEV222) dis- 
tribution function, based on Eq. (9) from [8] is 
shown in Eq. (10): 



LL{x, y, z, «i, au Pu "3, ":, ft, "3, »j, ft, m„ m^u /"h:) =/i | -rti In n^ 

^?[-('-(^)*)'*^'"('-(^)*h-']l 

+'!{S[-(lni* + ln'<,)+ln(l-(^'^^)fl.)"""*"' + ln(l-(^'^^)A)"^"'"' 



^'"(('-(^)fth*+(-(V)'^H 



ilnitfi - 2 



-'"(K-')^(('-(^)4**-('-(V)^)H ) 

-((.-(^)ft)-*.(i-(^)ft)-«n) 
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+ /3{-n3(lna, + Ina,+lna3) + |:[ln(l-(^^)j8. )"■'"'■-' 
+ln(l -(^)^.)""*- +In(l - (^)ft)'"'"^- 

+/4{i[-(lna, + lnai)+ln(l-(i=^)/3,)""'*-' + ln(l-(^')ft)'"«"^-' 
+ ln((l - (^)^,)'"'^'^+ (l - (^) A)"'«"^)''"'"" + ln((/nM- 1) 

((-(^)aH-(-(^)*M"')-(('-(^)'^H 



mt trivariatc association parameter 

nth], mw bivariate association parameter before and after the common period, respectively. 



HI «S 



2 i^f(Ph 8i) and 2 In/(r/, §5) take the form: 

1-1 i-l 



I [-lna,-(l-(^)/9.)''ft + ln(l-(^)A)"*-']. (II) 
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Similarly, ^ ln/(p/, ^/, Qi) and S ln/(rj, j,, gi) take 

the following form (bivariate relationship with both 
GEV marginals): 






mfc/A-l 



Vntb-l 



('-C^)«-H 



+'"(('"-0+(('-('^)ftH 



+ 



-(('-(^)«h* 



Given the complexity of the mathematical expres- 
sions in Eq. (10) and their partial derivatives with 
respect to the parameters* the constrained multi- 
variable Rosenbrock method, Kuester and Mize 
[10], was applied to obtain the maximum likelihood 
estimators for the parameters by the direct maxi- 
mization of Eq. (10). The required initial values of 
the parameters to start the optimization of Eq. (10) 
were provided by the univariate maximum likeli- 
hood estimators of the parameters for the case of 
the location, scale, and shape parameters. The Ini- 
tial values of the association parameters, bivariate 
and irivariatCt were set equal to 2, following the 
procedure developed by Escalante [8]. 



4. Reliability of Estimated Parameters 

The indicator selected to measure the reliability 
of estimated parameters when using the trivariate 



distribution as compared with the univariate coun- 
terpart was the asymptotic relative information 
ratio. 

Table 1 shows a sample of relative information 
ratios obtained by using the following set of 
parameters; 

ui = 15.0, ai =2,0, Pi^- 0.20 

1/2 = 12.0, 02=1.2, ft=-0.15 

Ii3=10.0, 03 = 1.0, /33=-0.10 



Table 1, Asymplotic relative information ratios of Ihe parame- 
ters of the TEV222 distribution for the maxima n3«25; m,**2; 
tttta — 2 



Parameter n^ 



"I 



ft 



U3 



^ 



^3 







25 



50 






1.086S 


1.3695 


1.5055 


25 


1.4460 


1J942 


1.6823 


50 


1.6295 


L7201 




15 


1-7408 









l.OHl 


1.2274 


1.3256 


25 


1.2712 


1.3753 


1.4356 


50 


1.3941 


1.4555 




75 


L4662 









1.2405 


13555 


1.4041 


25 


1^864 


1.4382 


1.4671 


50 


1.4514 


I.4S06 




75 


L4S83 









1.QS76 


1.3722 


1.5094 


25 


1.0599 


1.2263 


1.3333 


50 


LQ500 


L1683 




75 


1.0450 









1.0135 


I.2I51 


1.3065 


25 


0.9835 


1.1046 


1.1788 


50 


0.9734 


1.0604 




75 


0.9684 









1.2442 


1.3469 


13901 


25 


1.1324 


1.2069 


1.2507 


50 


1.0934 


1.1507 




75 


1.0736 









1.0882 


K0462 


1.0302 


25 


1.0604 


1.0417 


1.0313 


50 


1.0506 


1,0390 




75 


J.0736 









1.0126 


0.9553 


0.9334 


25 


0.9822 


0.9549 


0.9395 


50 


0.9790 


0.9542 




75 


0.9669 









1.2437 


1.0268 


0.9535 


25 


I.13H 


1.0246 


0.9704 


50 


K0923 


1.0223 




75 


K0724 







75 



1.5856 



L3821 



1.4312 



1^903 



13587 



L4140 



1.0217 



0.9218 



0.9166 
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S. Case Study 

A region located in Northern Mexico, with a total 
of six gauging stations, was selected to apply the 
proposed methodology to the flood frequency anal- 
ysis. Tables 2-4 show the results of the application 
of the trivariate extreme value distributions for the 
maxima to the data recorded in such gauging 
stations. 

In order to compare the goodness of fit between 
the univariate and trivariate maximum likelihood 
estimates of the parameters in stations considered 
in the case study, the standard error of fit, as de- 
fined by Kite [H], was obtained and is displayed in 
Table 5, 



Table 2. Coix^lalion coefficients and relative sample size^ for 
the tripled of statioru for the case study 



Table 5. Standard errors of fit for gauging stations of case i^tudy 



Tripleis of 

stations 



Correlation Relative sample sizes 
coefficient n\ «2 /13 /t^ n^ 



Acatitan-Sta Cruz-lKpaiino 

Choix-Huitcs-Sn Francisca 



0,926 


9 


2 26 





[J 


0.&69 





14 18 


7 






Table 3* Univariate maximum Jlkclihotxl estimates of the 
parameters of the GEV distributions defined by the data of the 
gauging stations of the case study 



Station 


Location 


Scale 


Shape 


Acatitan 


576.21 


283.80 


-0,62 


Choix 


236.69 


130,15 


-0A2 


Huitcs 


1564,78 


97«,87 


-0.57 


Ixpalino 


772.57 


473.97 


-0.38 


Sn Francisco 


926.53 


53156 


^0,65 


Sta Cruz 


835.74 


440.23 


-0.40 



Table 4, Trivariate maximum likelihood estimates of the 
parameters of the TEV222 distribution defined by the data of 
the gauging stations of the case study 



Station 


Location 


Scale 


Shape 


Acatitan 


568.93 


269.44 


-0.64 


Choix 


220.85 


I2«.2y 


-ojy 


Huites 


160330 


1038.53 


-0.68 


Ixpalino 


795.03 


490,86 


-0.46 


Sn Francisco 


943.69 


540.73 


-0.67 


Sta Cruz 


850.97 


467.74 


^0.52 





Standard error 


of fll 




Univariiitc 


Trivariate 


Station 


(GEV) 


(TEV222) 


Acatitan 


244.40* 


253.90 


Choix 


87,70 


58.80* 


Huites 


1024.00 


S3L90' 


Ixpalino 


537.90 


393.00* 


Sn Francisco 


350.80" 


401.50 


Sta Cruz 


497.20 


259.60- 



' Minimum standard error of fiL 



6, Conclusions 

The logistic model for trivariate general extreme 
value distribution for the maxima has been 
proposed* Asymptotic and data base results suggest 
that the proposed model is a suitable option to be 
considered when performing flood frequency analy- 
sis. 
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